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The pintle injector used for a liquid rocket engine is a newly re-attracted injection system famous for its wide 
throttle ability with high efficiency. The pintle injector has many variations with complex inner structures due to 
its moving parts. In order to study the rotating flow near the injector tip, which was observed from the cold flow 
experiment using water and air, a numerical simulation was adopted and a verification of the numerical model 
was later conducted. For the verification process, three types of experimental data including velocity distributions 
of gas flows, spray angles and liquid distribution were all compared using simulated results. The numerical simu- 
lation was performed using a commercial simulation program with the Eulerian multiphase model and axisym- 
metric two dimensional grids. The maximum and minimum velocities of gas were within the acceptable range of 
agreement, however, the spray angles experienced up to 25% error when the momentum ratios were increased. 
The spray density distributions were quantitatively measured and had good agreement. As a result of this study, it 


was concluded that the simulation method was properly constructed to study specific flow characteristics of the 


pintle injector despite having the limitations of two dimensional and coarse grids. 


Keywords: Spray characteristics, Pintle injector, Simulation, Experiment, Liquid rocket engine 


Introduction 


The pintle injector is a newly re-attracted injection 
system for a liquid rocket engine possessing several ad- 
vantages such as availability of one injector unit instead 
of many injector elements and wide throttleability with 
high efficiency. These advantages make the pintle injec- 
tor highly suitable for deep throttling of a liquid rocket 
engine due to its good efficiency at a wide range of throt- 
tling levels [1-5]. Unlike traditional coaxial injectors, the 
pintle injector has a complex shape due to its moving 
parts. For this reason, the flow characteristics near the 
injector outlet are very different from that of the coaxial 
injector. Therefore, fundamental studies of the pintle in- 
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jector are required to analyze spray characteristics of 
each specific injector. Son et al. performed such a study 
on spray angles of the pintle injector using a moving pin- 
tle type injector [6], however, it was difficult to analyze 
the detailed breakup and mixing characteristics through 
experiment only. From Son’s experimental research, an 
unusual spray structure was also observed near the pintle 
tip at specific injection conditions. Yue et al. conducted 
numerical combustion simulations of the pintle injector 
to investigate further [7]. During the investigates, al- 
though unique flow fields were visualized, the analysis of 
the detailed spray characteristics was not enough to 
properly study the effects of the injector shapes on the 
spray structure. 
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Nomenclature 

D diameter (mm) 

dopen pintle opening distance (m) 
L 

M momentum ratio 

m mass flow rate (kg/s) 
t thickness (mm) 

V velocity (m/s) 

We Weber number 
Greek letters 

P density (kg/m?) 

o surface tension (N/m) 


To overcome this challenge, in this study, in order to 
analyze the microscopic flow characteristics near the tip 
region, a numerical method has been adopted and verifi- 
cation was also conducted by comparing the contour with 
spray images from the experiment. Two dimensional do- 
mains with coarse meshes were used only for the phe- 
nomenological study. 


Experimental Data for Verification 


Experimental Setup and Conditions 


The pintle injector has coaxial injection ports with a 
moving pintle at the center port as shown in Figure 1. 
The pintle consists of a straight rod and cone shaped pin- 
tle tip with a constant spreading angle. The annular gap 
has a constant injection area and this area was fixed in 
this study. Detail specifications are arranged in Table 1. 
The center post protrudes from the injector plate and the 
diameter of the pintle tip is larger than the diameter of the 
center post. For this reason, the injector has a ring-shaped 
space at the post end with a 5.0 mm length and 4.0 mm 
depth. The flow rate and spray structure was adjusted 
according to the opening distance between the pintle tip 
and the center post. Although real propellants should be 
used to study spray characteristics, simulants, such as 
water and air, were used instead of liquid oxidizer and 
gaseous fuel. Water was supplied from the piston pump 
via two bypass valves for flow rate control and air came 
from the compressor as shown in Figure 2. In order to 
properly visualize the spray, a scattering technique was 
applied. The camera was installed horizontally at 45 de- 
grees for the axis and was constructed by the stroboscope 
and spray. Using this setup, the time averaged images of 
sprays were obtained and spray angles were measured. In 
this study, the measured spray angles were used to verify 
the simulation results. 

Additionally, one more experimental data was col- 
lected form a spray patternator as shown in figure 3. This 
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Fig. 1 Schematic of pintle injector 


Outer body diameter (mm) 
Inner body diameter (mm) 
Center gap diameter (mm) 
Center post diameter (mm) 
Pintle rod diameter (mm) 
Pintle tip diameter (mm) 
Center post length (mm) 
Pintle tip thickness (mm) 


Pintle end thickness (mm) 


Table 1 Dimensions of the pintle injector 


Do 13.5 
Di 12.0 
Deg 4.0 
De 8.0 
Dp: 3.0 
Dy 12.0 
Lp: 5.0 
tp 3.0 
tpe 1.0 


Bypass valves 


Fig.2 Schematic of spray visualization system 
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patternator consisted of 7 by 7 cylinders which each had 
a 20mm diameter. Spray was injected vertically into the 
patternator at a 50mm axial distance. The liquid in each 
cylinder was measured using an ultrasonic distance sen- 
sor. Measured mass distributions were normalized by 
collected total mass. For the verification, the mass con- 
centrations of 4 cylinders from the center point to the 
outside were compared with similarly calculated mass 
concentration. 

Four experimental cases at different injection condi- 
tions were compared as shown in Table 2. Case 1 was the 
reference and the pressure drop of liquid injection in- 
creased for case 2. Case 3 had a longer distance than that 
of case 2 and case 4 was designed to study the effect of 
the gas pressure drop. In order to generalize the results, a 
momentum ratio and Weber number were used as 
non-dimensional numbers shown in equation (1) and (2). 
Only bulk velocities at the injector exit were used in the 
numbers and the opening distance was adopted for cha- 
racteristic length of Weber number. 


y 
M = M gas gas (1) 
m a tig 
P gas (Vine —Viig i dopen 
Wee e ee (2) 
Oliq 
Indicator 
Ultrasonic 
sensor 
oN 
= 
Fig. 3 Schematic of spray patterning system 
Table 2 Experimental cases 
Case 1 Case 2 Case 3 Case 4 
aes eae 0.20 0.20 060 0.60 
distance (mm) 
Pressure drop of 
liquid injection 0.10 0.50 0.50 0.50 
(barg) 
Pressure drop of 
gas injection 0.05 0.05 0.05 0.20 
(barg) 
Momentum ratio 2.02 0.79 0.77 2.68 
Weber number 10.62 9.35 32.19 121.30 
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Experimental Results 


The scattering spray images of the four cases produced 
from the experiment are illustrated in Figure 4. At similar 
Weber numbers, the spray angle of case 2 was larger than 
that of the angle of case 1 due to an increase of momen- 
tum ratio. In addition, at the similar momentum ratio 
condition, although the overall spray structures of cases 2 
and 3 were similar, an additional spray layer was ob- 
served in case 3. Finally, the boundary of the spray layer 
was smeared in case 4. Only for case 4 was there a splash 
near the center post region and this phenomenon was 
consistently observed in the cases with high Weber num- 
bers. From these results, it could be deduced that the post 
shape strongly affected to the spray structure. In order to 
study this effect in detail, numerical method was adopted 
and verified. 


Results of Numerical Simulation 


Numerical Method and Conditions 


Numerical analysis was performed on the 2D axi- 
symmetric domain shown in Figure 5, using a commer- 
cial program. All boundary conditions were similar to the 
actual experimental conditions. The number of basic gr- 
ids was about 300,000 and refinements were conducted 
three times after convergence. An unsteady pressure 
based solver was used with the first order implicit me- 
thod for time and the second order discretization for 
space. The realizable k-epsilon model using standard 
wall function was applied to create a viscous flow model. 
The Eulerian multiphase model was adopted to simulate 
liquid and gas injection. The initial conditions and boun- 
dary conditions in all four cases are shown in Table 3 and 
these conditions corresponded to the experimental condi- 
tions. Injection temperatures were fixed to 300K for wa- 
ter and air in all four cases. At all outlet boundaries, at- 
mospheric conditions were assumed. 


Numerical Results 


In order to study the global spray structure, volume 
fractions of water phase were illustrated in log scales as 
shown in Figure 6. Similar to our experimental result, the 
sprays of cases 1 and 4 were strongly affected and mixed 
by gas flows while cases 2 and 3 were only lightly af- 
fected. There is a noticeable point for case 1 and case 4 
where the recirculation zones were constructed under the 
pintle tip. This rotating region could be the flame holder 
and stabilizer in combustion, however, the recirculation 
was unstable. And this recirculation will also cause a 
cooling problem in the pintle tip. Another remarkable 
point was the additional recirculation zone near the cen- 
ter post as shown in Figure 7. Streamlines colored by 
turbulent kinetic energies at the post region were differ- 
ent from those at the spray structures. Cases 2 and 3 had 
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(c) case 3 


(b) case 2 


(d) case 4 


Fig.4 Spray images from experiment 
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Fig. 5 Specification of numerical domain and boundary con- 
ditions 


two rotating flows while cases 1 and 4 had one large cir- 
culation. Strong rotating flows at cases 1 and 4 pulled up 
the liquid drops from the sheet and this phenomenon was 
also clearly observed in the experiment of case 4. Be- 
cause the recirculation at this region was more stable 
than that of the recirculation under the pintle tip, it could 
possibly be a flame anchoring point at the combustion 
process. 


Table 3 Numerical simulation conditions 


Liquid mass Gas mass flux Pintle opening 


flux (Kg/s.m°) (Kg/s.m°) distance (mm) 
Case 1 2255.45 77.92 0.2 
Case 2 3605.09 77.92 0.2 
Case 3 6195.22 77.92 0.6 
Case 4 6195.22 150.20 0.6 


Verification of Numerical Method 


For the purpose of verification of gas flow, two cases 
with only gas injection were compared with the results of 
the simulations as shown in Figure 8. The figure shows 
velocities along the radial direction at the axial point of 
50mm distance. The initial and final velocities from the 
experiment were very similar to the simulation results but 
the calculated profiles were more slender than the expe- 
riment results. Considering the limitations of the two 
dimensional simulation, these results were enough to 
analyze the phenomenological structures. 

From the contour images, numerical spray angles were 
extracted and compared with experimental results as 
shown in Table 4. Comprehensively, it was similar to that 
of the maximum spray angle as shown in case 3 and the 
minimum angle in case 1. This means that the overall 
spray angles have the same trends as found in the expe- 
rimental results. However, there was about a 25% error 
rate in the spray angles for the cases with high momen 
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Fig. 7 Contours of turbulent kinetic energy near the post region 
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tum ratios such as those in case 1 and 4. In actuality, the 
measurements were too obscured in order to properly 
define spray boundaries and additionally, the numerical 
grids were too coarse to capture the small droplets. For 
these reasons, a comparison of more precise quantities 
was required and therefore the data of spray density dis- 
tributions were compared. 

This method was more accurate than using that of the 
spray angle since the spray density distribution could be 
measured with exact quantity. While conducting the ex- 
periment, the distributions were measured using 49 cy- 
linders. For the axisymmetric assumption, measured liq- 
uid volumes of 3 cylinders and a single half center cy- 
linder were selected as shown in Figure 9. 

In order to apply the same collecting condition to the 
simulation, volume fractions at the same axial point were 
measured and grid data was integrated into 4 points as a 
single unit cylinder of 20 mm diameter. In cases 1 and 4, 
because of the high momentum ratios, the sprays were 
concentrated towards the center. The overall agreement 
of results was acceptable except for case 3 which was due 
to the limitation of breakup modeling in the 2-D domain. 
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Fig. 8 Gas velocities along the radial positions 


Table 4 Comparison of spray half angles 


Angle (°) Error 
Calculation Experiment (%) 
Case 1 21.1 28.4 25.7 
Case 2 44.6 47.2 5.5 
Case 3 53.8 58.1 7.4 
Case 4 39.5 52.4 24.6 


Fig.9 Selected patternator region in experiment 
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Fig. 10 Spray density distribution at radial direction 
Conclusions 


In order to verify the simulation method used for the 
spray characteristics study of a pintle injector, the results 
from our experiment and actual simulation were com- 
pared. Using three different sets of data, verification were 
performed on; (1) velocity distributions of gas flows, (2) 
spray angles and (3) liquid distribution. From these 
comparisons, the results can be concluded as follows; 

(1) Comparing gas velocities along a radial direction, 
the maximum and minimum values were well matched. 
However, the simulated velocity profile was more slender 
than the measured profile. 

(2) Spray angles of less atomized cases were finely 
modeled and recirculation structures were revealed. 
However, the highly atomized cases had an error rate of 
about 25% due to difficulty of measurement from the 
spray images. 

(3) Finally, a more precise quantity and spray density 
distribution were compared and the result showed rea- 
sonable agreement. 
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